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00 , Abstract. We study vortex nucleation for minimizers of a Ginzburg-Landau 

^ _ ^ ■ energy with discontinuous constraint. For applied magnetic fields comparable 

^ ^ ' with the first critical field of vortex nucleation, we determine the limiting 

Cn [ vorticities. 

1. Introduction and main results 

OO . In the framework of the Ginzburg-Landau theory, it is proposed to model the 

energy of an inhomogeneous superconducting sample by means of the following 
p I ■ functional (see [3llT0J) 

^. : (1.1) g,,H{i',A)^J^^(^\iV-iA)ij\' + ^(pix)-\if)' + \cmlA-Hf^ dx. 

^2 ' Here, fi C M^ is the 2-D cross section of the superconducting sample, assumed to 

occupy a cylinder of infinite height. The complex- valued function tp £ H^{Q;C) 
is called the 'order parameter', whose modulus jf/'P measures the density of the 
superconducting electron Cooper pairs (hence ip = corresponds to a normal state), 
^ ■ and the real vector field A = (^1,^2) is called the 'magnetic potential', such that 

C^~« I the induced magnetic field in the sample corresponds to curl A. 

C^^ ■ The functional (|l.ip depends on many parameters: ^ = k is a characteristic of 

the superconducting sample (a temperature independent quantity), H > is the 
intensity of the applied magnetic field (assumed constant and parallel to the axis 

C^ I of the superconducting sample), p{x) is a positive function modeling the impurities 

in the sample, whose values are temperature dependent. The positive sign of the 

^\ I function p means that the temperature remains below the critical temperature of 

the superconducting sample. 
In this paper, the function p is a step function. We take 5*1 an open smooth set 

S^ i such that 5*1 C r2, 5*2 = ri \ 5*1, and 



where a G IR+ \ {1} is a given constant. 

The above choice of p has two physical interpretations (see [6tll0j): 

• Si and S2 correspond to two superconducting samples with different critical 
temperatures; 

• The superconducting sample fl is subject to two different temperatures in 
the regions Si and 5*2, which may happen by cold or heat working S'2. 

Minimization of the functional l|l.ip will take place in the space 
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That is we do not assume a priori boundary conditions for admissible configura- 
tions, but minimizers satisfy natural boundary conditions. We study nucleation of 
vortices as the apphed magnetic field varies in such a manner that, 

lim = A , A e IR+ , 

e^o |lne| 

and we obtain that their behavior is strongly dependent on the parameter a, see 
Theorems 11.11 and 11.21 below. 

The singularity of the potential p causes a singularity in the energy as e ^- of the 
order of e~^, whereas the energy due to the presence of n vortices typically diverges 
like nine. In order to separate the different singularities of the energy, we let u^ be 
the (unique) positive minimizer of l|l.ip when H = (see Theorem l2.ip . the energy 
of Ue being of the order e~^. Then, if (-0, A) is a minimizing configuration of l|l.ip . 
it holds that (see Lemma [273)) . 

GeMii^, H) = GefiiUe^O) + Te,H { — .A 



and the configuration I — ,A\ minimizes the functional T. h introduced below, 

(1.3) :F,^H{^,A) = l^^(^ul\{S/-iA)^\' + ^{l~\^\'r + \cm\A-H\'^^ dx . 



The leading order term of the minimizing energy of l|1.3p will be described by means 
of an auxiliary energy introduced in l|1.5p below. Let A4{n) be the space of bounded 
Radon measures on fl, i.e. the topological dual of Cq (ri). A measure ^ S ^A{fl) can 
be represented canonically as a difference of two positive measures, ^i ~ p+ — pL-. 
The total variation and the norm of p, denoted respectively by |/i| and l|/il|, are by 
definition |/i| = /i_|. + /i_ and |j^|| = Ia^K^)- 

Given A > 0, we introduce an energy E\ defined on M{Vl.) n H^'^{Q) as follows. 
For pL£M{n)r^ H-~^{n), let h^, e H'^{n) be the solution of 

(14) I -div (^V/i^) + /i^ = Ai inf7, 

[ hf^ = 1 on dil, , 
where p is introduced in l|1.2p . Now E\{p) is by definition 



/ip-ll' 1 dx. 



(1.5) E^{fi)^jJj{x)\p\dx + J^^(^-^\Vh,\' + \ 

Theorem 1.1. Given A > 0, assume that 

TT 

Hm = A , 

£^0 |lne| 

then 

in the sense of T -convergence. Here the functional J- e,H is defined in lll.3\} above. 



The convergence in Theorem 11.11 is precisely described in Propositions 13.11 and 
14.11 below. 

Minimizers of l|1.5p can be characterized by means of minimizers of the following 
problem, 
(1.6) 

min I (^ -div(—-\/h] +h+l . ,. 

heH'o(^^) Jn\ >^ \P[x) ) p[x) 

-div(^v/i)+heA4(o) 



i™i%i<.ndx. 
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The above functional being strictly convex and continuous, it admits a unique 
minimizer, and so the functional Ex- 

Therefore, as a corollary of Theorem 11.11 we may describe the limiting vorticity 
measure in terms of the minimizer of the limiting energy E\ . 

Theorem 1.2. Under the hypothesis of Theorem ] l.R if {ip^^A^) is a minimizer of 
) [J.3|) . then, denoting by 

hi; =^ curlAj , fi^ ^ h^ + curl(i(^e , (V — iAi;)ipi,) , 

the 'induced magnetic field' and 'vorticity measure' respectively, the following con- 
vergence hold, 

(1.7) ^^^* in (CO'^(r!))* for all 7 e (0, 1) , 

(1.8) — -^ hf,, weakly in Hl{n) and strongly in W^'P{n) , V p < 2 . 
H 

Here ^^, is the unique minimizer of Ex. 

Sketch of Proof. 

Let us briefly describe the main points of the proofs of Theorems 11.11 and 11.21 ^ind 

thus explain what stands behind their statements. 

The starting point is the analysis of minimizers of l|l.ip when H = 0. In this case, 

(|l.ip has, up to a gauge transformation, a unique minimizer (ue,0) where u^ is a 

positive real-valued function. The asymptotic decay of Ue as e -^ is obtained in 

Lemma 12.21 

When H > 0, let {ip,A) be a minimizer of l|l.ip . Inspired by Lassoued-Mironescu 

[9], we introduce a normalized densitVQ 

Then |(/?| < 1 and we are led to the analysis of the functional .?>,// ((^9, ^). The 
main difficulty stems from the boundary layer behavior of the weight u^ , which has 
rapid oscillations comparable to e~^ near the boundary of Si. The challenge is 
then to construct a test conflguration with vortices localized near dSi and having 
the right amount of energy. Technically, one does that via a Green's potential Ge 
which we were not able to give a good control of it, see Lemma [3?2l To overcome 
this difficulty, we construct vortices situated at a distance ^j^" ! away from dSi. 

Those vortices being expected to lie on a union 7 of closed curves, a nice L^ bound 
can be shown to hold for Ge{x,x) on 7. Then, using a result of [6], recalled in 
Lemma [33 this L^ bound provides us with a family of well separated points serv- 
ing as the centers of vortices for the test conflgurations, see Proposition 13.51 We 
use then suitable methods to estimate from above the energy of the constructed 
configuration, yielding in Proposition 13. II the first part of Theorem ll.il This rough 
analysis will be carried out through Section [31 

In Section lU we establish a lower bound of the energy, see Proposition 14.11 This 
yields the second part of Theorem 11.11 

The proof of Theorem 11.21 is a by-product of the proofs of Propositions 13.11 and 
\A.1\ see Remark [T3l Finally, we close the paper with a discussion concerning the 
minimization of the energy l|1.5p . 
Remarks on the notation. 

• The letters C,C,M, etc. will denote positive constants independent of e. 

• For n S N and X C K", \X\ denotes the Lebesgue measure of X. B[x,r) 
denotes the open ball in R" of radius r and center x. 



Notice that </3 and ip have the same vortices. 
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• (•, •) denotes the scalar product in C when identified with R^. 

• ^e.ni'-P^ ^i U) means that the energy density of [(p, A) is integrated only on 
C/C f7. 

• For a e (0,1) sufficiently small, 5f = {x G 5*1 : dist(a;, 951) > a} and 
S^ ^{xeS2 ■■ dist(x, dSi) > a}. 

• For two positive functions a{e) and 6(e), we write a{e) ^ 6(e) as e -^ to 

mean that lim -— — = 0. 

e— 6(e) 

2. Preliminary analysis of minimizers 

2.1. The case without applied magnetic field. This section is devoted to an 
analysis for minimizers of (|l.ip when the applied magnetic field H = 0. We follow 
closely similar results obtained in [2- 

We keep the notation introduced in Section 1. Upon taking A = and if = in 
(|l.ip . one is led to introduce the functional 

(2.1) g^^u):=l^^(^\yu\' + ^{p{x)-u'r'^ dx , 

defined for functions in H^{fl; M). 
We introduce 

(2.2) Co(e)= inf g,{u) . 

«effi(f2;K) 

The next theorem is an analogue of Theorem 1.1 in [7]. 

Theorem 2.1. Given a G R+ \ {1}, there exists Eq such that for all e g]0, eo[, the 
functional f2.1\) admits in H^{fl;W) a minimizer u^ S C^(S'i) U C^(S'2) such that 

min(l, Va) < w^ < max(l, \/a) in fl. 

If H = 0, minimizers of fl-l]) are gauge equivalent to the state (ue,0). 

We state also some estimates, taken from [71 Proposition 5.1], that describe the 
decay of Ug away from the boundary of Si . 

Lemma 2.2. Let fc G N. There exist positive constants Eq, S and C such that, for 

all e e]0, eo], 
(2.3) 

C 



,, , , 6dist{x,dSi) 

( 1 — Ug) exp ' 



miSi 



, r- X / (5dist (a;, 95*1)1 

\\J a — Ug) exp ' 



ff'=(S'2) ^ 



Finally, we mention without proof that the energy Co(e) (cf. (|2.2p ) has the order 
of e~^, and we refer to the methods in (Tj Section 6] which permit to obtain the 
leading order asymptotic expansion 

Co(e) = £iM + c2(a)+o(l), (e -> 0), 
e 

where ci(a) and 02(0) are positive explicit constants. 

2.2. The case with magnetic field. This section is devoted to a preliminary 
analysis of the minimizers of (|l.ip when H y^ 0. The main point that we shall show 
is how to extract the singular term Co(e) (cf. I|2.2p ) from the energy of a minimizer. 
Notice that the existence of minimizers is standard starting from a minimizing 
sequence (cf. e.g. [4]). A standard choice of gauge permits one to assume that the 
magnetic potential satisfies 

(2.4) divA = infi, n ■ A = Q on 9^^, 
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where n is the outward unit normal vector of dfl. 

With this choice of gauge, one is able to prove (since the boundaries of il and 5*1 
are smooth) that a minimizer (^, A) is in C^{fl;C) x C^(fl;M.'^). One has also the 
following regularity (cf. [3 Appendix A]), 

V' e C^(Si;C) U C^(S2;C), A e C^(Si;m.^) U C^(S2;M.^). 

The next lemma is inspired from the work of Lassoued-Mironescu (cf. [9]). 

Lemma 2.3. Let {4',A) be a minimizer of ^.1]) . Then < {tpl < u^ in Vt, where 

Me is the positive minimizer of h2.1\) . 

Moreover, putting ip ~ ^, then the energy functional fl.l]) splits in the form : 

(2.5) g,.H{i^,A) = Co{e)+T,,H{v,A), 

where Co{e) has been introduced in f2.2\} and the new functional Te,H is defined 
by : 

T.^hW^A) = ^ Ll\{^ - zA)(^p + -^u\{\ - l^n^ + Icurl A - H\^ Ax. 

3. Upper bound of the energy 

3.1. Main result. The objective of this section is to establish the following upper 
bound. 

rr 

Proposition 3.1. Assume that lim = A with A > 0. Given u G MiVl) C\ 

e^o |lne| 

H~^{Q), there exists a family of configurations {{^Pe, ^e)}e it^- H^{^', C) x H^{fl; R^) 

such that \\ipe\\L=^(n) < Ij 

(3.1) ^^^^^^M in {Cf-my V7e(0,l), 

and 

\imsnpF,,H{(Pe,A,)<- [ p{xMdx+ [ ( _-|V/i^p + |/i^ - Ip ] dx . 

3.2. Analysis of a Green's potential. This section is devoted to an analy- 
sis of a Green's kernel, i.e. a fundamental solution of the differential operator 

— div I V 1+1. The existence and the properties of this function, taken from 

[UITS], are given in the next lemma. 



Lemma 3.2. For every y <E O. and e e]0,l], there exists a unique symmetric 
function VI'k Q.3 [x,y) ^^ Ge{x, y) G M_|- such that : 

(3 2) I -div (^-^V,G,{x,y)^+G,{x,y) = <S,(.t) in V'{n), 

I G,(x,y)|^^g,^-0. 

Moreover, G^ satisfies the following properties: 

(1) There exists a constant C > such that 

Ge(x,y)<C(|ln|x-y|| + l) , V (x, y) G H x II\ A , V e e]0, 1] , 

where A denotes the diagonal in M? . 

u (x) 

(2) The function Vg : fl x Q. 3 {x , y) i-^ Ge{x, y) ^ — ^ — In jx — y| is in the class 

G^{VL X rj;R). 

(3) Given q E [1, 2[, there exists a constant C > such that 

\\ve{-.v)\\w^.,[n)<C, yyeU, Vee]0,l]. 
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One feature of the function G^ is that we know its homogenized Umit Go- 

Lemma 3.3. Let 1 < g < 2. The function G^ converges weakly in W'^''^{Vl x Q) to 
Go the solution of 



-div ( zjzx^xGoix,y) ] + Go{x,y) = Sy{x) in I?'(r!), 



(3.3) <! \p{x) 

Moreover, G^ converges locally uniformly to Go inflxQ,\A where A is the 

Proof. Since Gg and Go are symmetric, it is sufficient to estabUsh weak convergence 

of Ge{-, y) for y G fl arbitrarly fixed. 

Lemma [3^21 assures that Gj(-,y) is bounded in W^''^{fl) uniformly with respect to 

y and e. Hence, Ge{-,y) converges weakly to a function Go(-,j/) in W^''^{fl). Let 

us prove that Go solves p.3p . 

Let f e G^{n). It is sufficient to prove that as e ^ 0, 

/ „, . \7xGe(x,y) ■ \/x'p(x)dx ^ / , . \7xGo(x,y)-\/x'fi(x)dx 
Geix,y)(p{x)dx ^ / Go{x,y)ip{x)dx. 



The second convergence is actually immediate since VK^'^(17) comapctly embedds 
in L'^{fl) and Ge{-,y) converges weakly to Go{-,y) in W^''^(fl). In order to prove 
the first convergence we notice that, 

1 /■ 1 

\/xGe{x,y) ■\/xf{x)dx ~ / ——-\/xGo{x,y)-'^x'f{x)dx 



nui[x) JnP{x) 

+ / —r-^^x{Qe{x,y)~Go{x,y))-^xV{x)dx, 

so it suffices to prove that 

/ I n. s ---]VxGeix,y)-Vx^ix)dx -^0 as e ^ . 

Using the uniform bound on Ge{-,y) in W^''^, we write, 
1 1 



n 



Ug(x) p{x) 



VxGeix,y) ■ VxVix) 



dx 



I 1 ^ ^ \^/P / r \ i/« 



1 1 \'._ _ . a"' 1 1 



Now, we use the asymptotic behavior of u^, which converges pointwise to p{x) and 
remains uniformly bounded from above and below. Therefore, Lebesgue conver- 
gence theorem yields the desired convergence. 

To complete the proof, we refer to [H Lemma 3.1] to find the uniform convergence 
of Gj to Go. D 

Next we cite a result from [HI Lemma 5.5]. 
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Lemma 3.4. Let {fe)E£]o.i] C C([0, 1],K+) be a family of continuous functions. 
Assume that there exists a constant C > such that 

\\fe\\LH[OS])<Clli\llie\, V£e]0,l]. 

There exist constants K > and cq g]0, 1[ such that, given a family {N{e)) C N 
satisfying N{e) ^ 1, there exists a family ((5(e)) c]0, 1[ and a sequence (im)meN C 
]0, 1[ and 

\fM)\ < K\n\\ns\ , ^ < |tf+i -tt\< S{e) + ^ , 

Vie{l,2,--- ,iV(e)}, Vee]0,l]. 

The next proposition provides us with points enjoying useful properties. These 
points will serve to be the centers of the vortices of the test configuration that we 
shall construct in the next section. 

Proposition 3.5. Let /i G CQ{fl) be continuous and compactly supported in fl. 
Assume that (fi(e))eG(o,i] C N is a family of integers such that ci|lne| < n{e) < 
C2e~^ for constants ci and C2 independent of e <E (0, 1] . 

// the restriction ji of ji to Si (respectively S2) is different from 0, then there exist 
constants c > 0, C > 0, a family of integers N{e) e N, points af in Si (respectively 
S2) and degrees df = ±1 such that: 

(1) N{e) = n{e){l + o(l)) as e ^ ; 

(2) \cq ~ a'^ > -^= > ce for all i^j; 

(3) dist (af , dSi) > ^^ "!"' and dist(af , 917) > c for all i ; 

(4) \v,{al,al)\ < |lne| for alii ; 

(5) -T-TT / Mi ^ 27r7-:— 7 in the weak sense of measures, for u5 any measure 

^ ' N{£) -^ ' ll^ll ' J M-i y 



supported in B{af,ce), of constant sign, and such that /if(rj) = 2TTdf. 

Proof. For i = 1,2, we denote by S'f = {x e S, : dist(x, dSi) > ^^^j\^}. We only 

present the construction in Sf as it holds exactly in S"!- 

We partition Sf into squares K of side-length ^(e) where £{e) is chosen such that, 

< e{e) < 1 ase^O. 
V'^(e) 

We denote by /C(e) the family of such squares that he entirely inside Sf. 
For K & IC{e), we set 

where 

M,= ^ \^iK)\. 

K<£K(e) 

Since /j, is continuous with compact support in Q., and (.{e) tends to as e ^ 0, we 

get that Me -^ \\jl\\ as e ^ 0. 

Now >, ^K = n, hence Lemma 7.4 of [11] provides us with nonnegative integers 

K£K(e) 

m^{K) such that 

(3.4) Y. ^MK) = nis), 

KeK{e) 
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and 

HK)\ 



(3.5) 



m^{K) - n- 



M, 



< 1. 



Since /i is bounded, we may find M > such that \^-i-\{K) < Ml{e)'^ for all K e K.{e). 
Consequently, we infer from l|3.5p that 

me{K) < 1 + nMt^ = OinMt^) , 

(recall that i<£<las£^0). 

Now we decompose /C(e) = /Ci(e) U /C2(£), where K £ /C2(e) if and only if 

dist{dK,dSi) > -^ where i? > is picked small enough so that each connected 

component of Si contains a disc of radius R. 

Now assume that K e /Ci(e). Let PeiK) be the least integer less than y/m^iK). 

Then we may pick Pe{K) points 6f lying on one side C oi K and such that 

(3.6) \bt-^\>^^> ^ 



Next, we can construct Pe{K) segments Lf in K, each passing by bf and perpen- 
dicular to the side C (hence each Li has length i{e)). 

Since K G /Ci(e), we may find closed smooth curves ji D Li in ft such that each 
-fi is the boundary of a simply connected set Ui <Z Vl and 2ttR < \ji\ < \dSi\. By 
this way, referring to [H Proof of (5.4)], we may get a constant C" > independent 
from e and such that, for all i, 



'7i 

In particular, it holds that 



/ |we(a;,a;)|dx<C"ln|ln£|. 



f \v,{x,x)\dx<C'\n\\ne\ 

JLi 



'Li 

Let us parameterize Li by (0,^(e)] 3 s i-^ x{s) e Li, and let us define the rescaled 
continuous function 

(0, l]3t^ f,{t) = V, {x{i{e)t),x{i{e)t)) , 

so that ||/e||Li(o,i) < C"(£(£))^^ ln| lne|. Invoking the result of Lemma [331 we 
may pick Pe{K) points (&f ,) C Li such that|3 

and 



Denoting by {a\) the family of the constructed points in all the segments Li, we 
have actually constructed [pe(X)]^ points (a|) C K such that 

1 , ^ ^M ^^;ln|lne| , . ., C 

\v,{al,al)\<K'^—l, l4-a|> 



Recalling the assumptions on n[e) and £(e), we get actually the two desired prop- 
erties (2) and (3) stated in Proposition 13.51 above, with the constant c > chosen 
sufficiently small that c < dist(supp/j, , dO). 

It remains now to continue the construction of the points (af ) filling the squares 
K e K,2[£)- Notice that Corollary 5.3 of [6] (more precisely an adjustment of its 



We apply the lemma with N{e) = \/n{e) so that we get ^n(e) > Pe{K) points. 
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proof) provides us with a constant C > such that, for all K G /C2(e) and e G (0, 1), 
we have, 

||ue(x,a;)||Lo=(i^) < C . 

Hence, it is sufficient to construct any well separated rae{K) points in this case. 
This is exactly the case of [HI p. 144]. Now, the integer iV(e) is defined as 

KelCiie) K<^K2(e) 



That N{e) — n(£)(l + o(l)) is due to the assumption we made on ^(e), \/n{e) ^ 

(.{e) < 1 as e ^ 0. 

Now, having constructed the family of points (af), we see that property (3) stated 

in Proposition [33] is just due to our construction of the points being in Si, and our 

choice of the constant c being so small that dist(supp/i , dVt) > c. 

Now we define the family (d^). If ^j,{K) > 0, we assign the degree di = 1 to each 

Qi G K, otherwise we associate the degree di = —1. The proof of the last property 

(5) in Proposition 13.51 is exactly as that given in [HI p. 145]. D 

Proposition 13.51 will be used in the following context. Let /i G Co(r2) be contin- 
uous and compactly supported in fi. Take positive integers ni(e) and 712(2) such 
that 

ci| lne| < ni{e) + n2[e) < C2e^^ , 

for positive constants ci and C2 independent from e. 

If the restrictions jl and /i of /z to Si and S2 respectively are both different from 0, 
then we get Ni{e) ^ ni{e) points (af) in Si and N2{e) ^ n2{e) points (bf) in S2 
satisfying properties (l)-(5) of Proposition [3751 In particular, we set 

af ViG{l,2,...,ni(e)}, 

bf \/ie{?ii{e) + l,n2{e) + 2,--- ,ni{e) + n2{e)} , 

and we define the measures //f by 

if a; G Bi^xfjCe) 



(3.7) ^iKx) 



,2 £.2 



c^e 



otherwise . 



If /i (respectively fl) is zero, we may still choose the points af (respectively bf) 
arbitrarily so that properties (l)-(4) of Proposition 13.51 are valid, and we define the 
corresponding measures fj,f to be zero by pure convention. 
With these notations, we get as an immediate consequence that 

Ni{6)+N2ie) 

(3.8) /^e := ^^ E 

n{e) ^-—' 

where 

(3.9) n{e) = 

Ni{e) + N2{e) otherwise. 

Lemma 3.6. Under the hypotheses and notations above, one has 
(3.10) 

limsup — r-) V/ Geix,y)dfil{x)dfi'{y) < -—r^ GQ{x,y)dn{x)dfi{y) . 

e^o \n(e)J f^.JB^xB, \\^^\r Jnxn 

Here, for all i, Bi denotes the ball B{xf,ce). 



L 


-'^M ^-(^o''^(")r' 


Niie) 


if supp yU C S*! , 


N2ie) 


if supp n d S2, 
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Proof. Given a > 0, let Aq, = {(a:,y), |x — y| < a}. We have the following 
decomposition, 

(3.11) 4l] / G,Ay.lAy.]= f G,d^led^ie + ^y2 f G,d^l|d^l'.. 

Knowing from Lemma [3.31 that G^ converges uniformly to Go in 17 x ri\AQ, we 
write using l|3.8p . 

47r2 



(3.12) 



lim 



Ged/igd/ig 



God/i d/i . 



/axn\A„ IIa^P Jnxn\A^ 

Let us estimate now the last term on the right hand side of p. lip . Since the 
supports of ^f and /i^ are disjoint for i ^ j, we may write by Lemma [3?2t 



„2 Z_^ 



i=^3 



G^d[i\d[fA 



Ao 



< — V 



n^ 



ln|a;-2/|||At-|(a;)|Mf|(y)dxd2/, 



where G > is a constant independent from a. and e. 



Proposition l3 . SI provides us that the points a\ are well separated, i.e. |af — a^| > 
This actually permits us to write (see [HI p. 147] for details). 



G 



(3.13) 



\Y.I Ged/.;MM5 < G / (I In |x - y| I + 1) dx dy . 

^ ■_, ■ J A„ J An, 



Substituting l(3TT2l) and fSTTSJ) in fSTTTj) . we deduce that, 

G,d^ll{x)d^l){y) < 

iBiXBi 



lim sup ^V" / Ged^l|{x)d^l){y) < ^^ / God^d^ 

£-0 n^f^^B^xB, hA\ Jnxn\A^ 



+G I (I ln|a;-2;|| + l)dxdy. 
Ja„ 



with a being arbitrary in the interval (0, 1). Making a ^ (recall that In \x ~ y\ is 
in L^), we obtain the desired bound of the lemma. D 

Definition 3.7. A family of points (a^) C fl satisfying Properties (2)-(4) stated in 
Proposition 1 3. 5\ is said to be a well- distributed family. 



Now, let /le : fl — >]0, 1[ be the solution of the equation: 



1 



(3.14) - div — V/i^ 1 + /i£ = inQ, h, ^ I on dD, , 
where u^ is introduced in Theorem [2TTJ We define also, 

(3.15) Me) = J (^|V/ieP + \h, - 1\A dx. 

Next, we state a remarkable energy-splitting due to Bethuel- Riviere [2]. We find 
it in [6l Lemma 5.7]. 

Lemma 3.8. Consider {ip,A) e H'^{n;C) x H'^{n;M.'^) and define 

A' = A-^jV^h,. 
uj. 

Then we have the decomposition of the energy, 

(3.16) T,.H{ip,A)^H^Me) + C,.H{^,A')+no + 2H /" (/i, - l)^((/p, A'), 
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where 

(3.17) no = H^ f \i\ip\''-l)\VK\^, fi{ip,A') = h' + cnrl{i^,{V-tA')^). 

Here, the functional C^^h is 

(3.18) C,.h{v, A') = ^ ul\{\/ zA')^P + IcurMf + ^4(1 - l^^f . 

Lemma 3.9. Given e G (0, 1), a well-distributed famili^ of n = ni + n2 points 
(ai) C ft (rii points in SI) together with degrees (di) C { — 1,1}, there exists a 
configuration iip,A) such that, /if being the uniform measure on dBi = dB{ai,ce) 
of mass 2Trdi, and letting 

n 

IJ'' =^l^t, 

i=l 

we have for some a. G [0, 1[, 

n 

Ce.ni'P, A') < ce"n| Inep + c'n o{\ lne|) + 2tt\ lne| ^p(a,) 

(3.19) *"^ 

T,,H (^, A) ^ C,.h{<p, A') + H^Me) + 2H f [h, - l)d^lix) 

(3.20) ^ 



+ O I neH + ne^H^ + {ri^eH + eH''-)^! C^^h{v. A'\ 
Moreover, for any Q < ^ <1, it holds that. 



(3.21) -M^,A) - m||co..(o) < CeWl + ei/ + ■ /•^^'«('^'^') 



Proof. We construct a test configuration {if, A). We define a function /i in fi by 
h = h' + Hh^ where h^ has been introduced in l|3.14p and h' is the solution of 

(3.22) I -div(^V/.')+/^' = M^ in f7, 

[ h' ^Q on dn. 

Note that h'{x) ~ J^ Geix,y)dfi^{y). As a consequence, we have 

(3.23) / f ^^ + \h'A dx^ f f G,ix,y)dfiix)df,iy). 
Jn V "e / Jn Jn 

Now we define an induced magnetic potential A = A' + ^V^/ie by taking simply 

curlA' = /i'. 

This choice is always possible as one can take A' = V^g with g G H^{Vt) such that 
Ag = h' . We turn now to define an order parameter Lp which we take in the form 

(3.24) ^ = pe*^ 
where p is defined by 

if X e {JiB[ai,ce), 

(3.25) P{x)^\ I '^ x(^\J,B{a,,2ce), 
— — ^-1 if 3i s.t. X £B{ai,2ce)\B{a,,c£). 



See Definition 13.71 above . 
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The phase 4> is defined (modulo 27r) by the relation: 

(3.26) V0-A' = -\I^h' in^\\JiB{a„ce), 

and we emphasize here that we do not need to define cj) in regions where p vanishes. 
Having defined {'^^A) as above, we estimate Ce^ni'^iA'). Recall that 

(3.27) C,.h{^, J^) = ^^ [ul\Wp\'' + p'ul\W^ - Af + \h'\^ + ^(1 - p^f^ dx. 
From (|3.25p and using the uniform upper bound of u^, it follows easily that, 

(3.28) ^ {nl\Wp? + ^(1 - P^f^ d.T < Cn. 
Thanks to l(335| . (p^ejl and the definition o{ h = h' + Hhe, we have. 

Replacing this in p.27p and invoking (|3.28p . we find, 

(3.29) C,Mv,A')<jO-^^^^^^^^^^^ + \h-HhA dx + 0{n). 
Thanks to l|3.23p . we deduce the upper bound, 

(3.30) C,,„{^,A')< [ G,dp{x)dp{y)+0{n). 

JOxO 

We now decompose the double integral in two, 
(3.31) 

/ G,dp'{x)dp^y)=J2 [ GApK^Wl{y)+Y, [ GApl{xW]{y). 

2 

Let us estimate the first term in the right hand side. Writing G^ = v^ — ^ln \x — y\, 

we have 

(3.32) 

G,ix,y)dptix)df,Ky) = f (v,ix,y) -^\n\x-y\) d^^f (x)d/.f (j/), 

Assuming {x,y) £ supp/i^ x supp/i^, we get 
ul{x) In \x - y\dp'i{x)Apl{y) 



BiXBi 



/ / u2(a, + c£e'''i)ln|cee'^i -C£e'^^|d6lid6l2 
Jo Jo 



2tt 



27rlne / ui{ai + cse"^')de 







Here and in the sequel, Cg or Cg denote constants bounded (uniformly in e) from 
below and above, and that may change from one line to another. 
The points a^ being away from the boundary of Si, Lemma 12.21 gives that 
Mg(ai + cee'^^) is exponentially close to p{ai). Thus, we may write, upon con- 
sidering the summation, 

(3.33) Y. ^^\n\x-yWl{x)dpl{y) = C,n + \neY,p{a,). 

, , JB,xBi ^T^ , n 
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Let US now estimate ^"^-^ / \v^{x,y)\diif{x)dfi^{y). Referring to Corollary 5.3 

JBiXBi 

Cq, 



iBiXB. 

in [6], we know that for some a G (0, 1) and all 77 G (0, 1), 



\\Vei-,y)\\cO'''{{xen : dist(a;,aSi)>j;} ^ ^2 



< 



Consequently, the following estimate holds, 
(3.34) 

J2 [ \yeix,y)\d^l|ix)d^lUy)<J2 [ (\v,ia,,a,)\+C^]d^i|ix)d^l|iy) 

„_, JBiXBi :_-, JBiXBi \ V / 



4=1 ■ 



47r Cn^r + An n\ve(ai,ai)\, 



.,, 21n|lne| 

With r] = -^w^^- 



By our hypotheses, we know that \ve{ai, ai)\ <^ | Inej. Thus, 

(3.35) V/ \ve {x, y)\dnUx)dfiUy) <ce°'n\ Inef +c'no{\\ne\). 

,_i JBiXBi 



i=l 



Combining ([S^S]) together with ([S^S]) . we get ([SlQ]) . 

The proof of the properties l|3.20p - l|3.2ip are exactly as that given in |1H p. 140- 

142]. D 

3.3. Proof of proposition 13.1], completed. We first assume ^ ^ is a continu- 
ous and compactly supported function. Let rii = [:||-|^|(5'i)] and 712 = [57ImI('5'2)] 
where [•] denotes the integer part. We take n = rii + n2. Since Si and 5*2 are 
disjoint and cover fl, we have H^tH = |/i|(ri) = |a«|('5'i) + |mI('5'2) and consequently. 



Proposition l3 . 51 provides us with a well-distributed family of points that serves as an 
input in Lemma [3^ Thus, we get configurations {(p^, A^) and associated measures 
jif such that 

I "i+"2 , . * 

^e := — Y^ a4 ^ M . weakly in [Cg'^j , 
i=i 
and the estimates IpTTO)) . ^J^-^^ and i(33T|) being all valid. 
Since H has the order of | Inej, Lemma [3^ and estimates p.l9p - p.36p yield. 



C,,H{^,A')<2n\\ne\J2pia^) + 'i^^7r^ / Goix,y)dfi{x)d^iiy) + oiH^), 
^ Wf^W JnJn 

Inserting the particular choices of ni, 712, n = rii + 712 and using the definition oi p 

being piecewise constant, we get, 

(3.37) 

C,.,Hiip,A')<H\lns\ f p{xM{x)dx + H^ f j Go{x,y)d^l{x)d^i{y) + o{H^). 

Again the hypothesis on H = 0(| lne|) gives that the remainder terms in p.20p are 
0(1), leading thus to 

^e,Hiv.^) < H^Joie) + H\lne\ f pix)\fi\{x)dx 

Jn 

+2H^ f {h, - l)d^i, + H^ I I Go{x, y)d^i{x)d^Iiy) + o{H^). 
Jn Jn Jn 
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Using l|3.15p . we rewrite the preceding formula in the following explicit form, 
(3.38) 

T,.h{v,A)<H^ f (\Ylhl + \h,~l\A dx + H\\ne\ j p{xM{x)Ax 
Jn V "e / Jn 

+ 2i/2 f (^h, - l)d^i + H^ f G,(x, y)M^)My) 

+ 2ij2 f (he - l)(d^e - d^l) +H^ f (Go - Ge)d^d^ + o{H^). 
Jn Jnxn 

Let us define the functions Ufj,,^ and h^,^ by 

(3.39) { -div(^VC/,,,)+C/,,, = A. inf7, 



C/^,e = on dn. 



(3.40) 



1 



-div ^-V/i^^e + /i^^e = fi in fi, 



/i^e = 1 on dfl. 



Remarking that C/p,e(a:) = / Ge{x,y)d^{y), we get as a consequence 

Jn 



(3.41) I l\YEl^ + \u,,^^\-]dx= I I G,ix,y)dfiix)dt,{y). 



P^el , \TJ |2 



Writing /i^,e — 1 = [/^.^ + /i^ — 1 and replacing l|3.4ip in in (|3.38p leads, after some 
calculations to 



+ 2 / (/le - l)(dMs -dfi)+ [ (Go - Ge)d^dAi + o(l). 
Jn Joxn 

We define ft.^ by 

(3.43) I -^^"(^^'^'^)+^-='^ ^""' 

[ /i^ = 1 on 9r2. 

By a standard compactness argument similar to that given to Lemma [33j we check 
that. 

Using a Cauchy-Schwarz inequality, it follows from Lemma [373l and the boundedness 
of /i that 

(3.45) lim f {Go-Ge)dfidn = 0. 

Next, noticing that H^ie — l||cO(n) ^ 1) it follows immediately from the convergence 
of Ate to^in (G^(r2))*, 

(3.46) lim / (/le - l)(dAie ~ dfi) = . 

=-~*o./n 

TT 

Replacing ((Oil) - ((06)) in ((3^ then using lim r = A, we get finally. 



-imsup ^--^^^'^^ < \ f p{xM{x)dx+ f f^\S/h,\' + \h, 

e^O H^ X Jf^ Ji^ \p{x) 
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Moreover, since l|3.2ip holds and L^Mi'^T^') — C\\ne\'^, we have 



1||M^, A) A-ll(cr(n)). < Ce-^ (^1 + V ^^^^^ j ^ °(1)- 

We conclude that l|3.ip holds, which finishes the proof in the case where /i is a 
continuous and compactly supported function. The general case where /x € M{fl)r\ 
H^^{n) follows by a standard approximation argument, see [TH p. 149] for details. 

4. Lower Bound 

4.1. Main result. The objective of this section is to prove the lower bound stated 

in Proposition 14.11 below. 

Given a family of configurations {{(pe^A^)}, we denote by 

(4.1) j,^{i(p^,{W~iA^)ipe), h,^ curl A,. 

TT 

Proposition 4.1. Assume that lim t- r = A with A > 0. Let {((pe,Ag)}n be a 

e^o I lne| 

family of configurations satisfying Te^HifcA,.) < CH^ and \\^Pe\\L'^(n) < 1 for a 

given constant C > 0. 

Then, up to the extraction of a subsequence e„ converging to 0, one has, 

V 7^(0,1) ^'^^^"/^"^ ^M in {C''\n))\ 

i h 

-jf-^J, -jf^^ weakly in L2(0) . 



Moreover, fj. ~ curl j + h and 

(4.2) 



lifJo^f :^^^ > E,U.) + J^ La:) 



J + ^V^/ip 
p{x) 



h^\'\ dx. 



Here, the energy E\ and the function /i^ are introduced in ( [J.5)) and {1-4^ respec- 
tively. 

4.2. Vortex-balls. In this section we construct suitable 'vortex-balls' providing a 
lower bound of the energy of minimizers of Ijl.ip . Recall the decomposition of the 
energy in Lemma [273l which permits us to work with the 'reduced energy functional' 

Notice that, by using (u^, 0) as a test configuration for the function l|l.ip . we deduce 
an upper bound of the form : 

(4.3) Te,H{^,A)<CH^, 

where ip = if^/u^, {ip, A) a minimizer of Ijl.ip . and C > a positive constant. 

We recall the hypothesis that there exists a positive constant C > such that the 

applied magnetic field H satisfies 

(4.4) i7<C|ln£|. 

The upper bound l|4.3p provides us, as in [IT], with the construction of suitable 
'vortex-balls'. 



Proposition 4.2. Assume the hypotheses i4.4\ )- Given an open set U G fl and a 
number p e]l, 2[, there exists a constant C > and a finite family of disjoint balls 
{B{{ai,ri)}i^j such that, {ip,A) being a configuration satisfying the bound i^TSp, 
the following properties hold: 

(1) B{ai, ri) C U for all i ; 

{2)w = {x£U : \^{x)\<l-\\ne\-^}c\jB{a„r,). 



16 
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(3) 5^r, <C|lne| 



-10 



iei 



(4) Letting di be the degree of the function yy/\(p\ restricted to dB{ai,ri) if 
B{ai,ri) C fl and di = otherwise, then we have: 

(4.5) / ul\{\/ - iA)ip\^ dx + f \cnylA-H\^dx> 

JB{ai,ri)\u} JB{ai,ri) 



2TT\d,\[ min u^ ) (|lne| -Cln|ln£| 

yB(a,.ri) 



(5) 



277 ^ diSai - curl(A + {iip, V a"^)) 



iei 



<C|lne| 



IVo' '"C^) 



We follow the usual terminology and call the balls constructed in Proposition [47 
'vortex-balls'. The proof of Proposition 14.21 is actually a simple consequence of the 
analysis of [TT] . 

4.3. Proof of Proposition 14. IL Let us consider smooth and open sets Ui C 5*1 
and U2 C S2, and let us denote by U their union. 

Applying Proposition 14. 21 in Ui and U2 respectively, we get families of balls B{ai, ri) 
and degrees di such that, 



Te,H{'^e.Ae,Ve) > 27rJi^ V 



m 



H 



min unM + o(i), 



B(ai,ri) 



where V^ ~ UiB{ai,ri). 

Recall that u^ converges uniformly to the function p in U. Thus, we may rewrite 

the above lower bound in the following form, 



TeMJ^e^A^V,) ^g^ llnel /^ Eg.gai Md , „ Ea.gl., M» 



iJ2 



H 



H 



H 



0(1), 



Using the bound ||^e||L~(n) < 1, we write |(V-a,)(^eP > |(Pe|'|(V-ae)|2 > |j,|2. 
Consequently, we have. 



(4.7) 



ff2 



> 



"? 



u\v. 



H 



dx . 



With U = fl, we infer from the bound J^^^h{Ve,As) < CH^ that up to the extrac- 
tion of a subsequence, -j- and — j respectively converge weakly to j and h in L'^{Vt). 

Thanks again to the bound (|4.3p , we deduce that — °' ,^^ and — °' ^f are 



H H 

bounded. Hence the measures — — — -j and — — — -^ are weakly compact 

in the sense of measures, and thus, up to extraction of subsequences, they respec- 
tively converge to measures ^1 and /i2 in {Cq{Ui))* and {Cq{U2))* ■ Thanks to the 
last property of Proposition 14.21 we get upon setting /i = curl j + h, 



A'l = M|„i ' M2 = M|„2 ■ 

Therefore, combining (|4.6p - (|4.7p and using the uniform convergence of u^ to p in 
[/, we deduce that, 

Te,H{^e,A„U)_ ^ Te,H{^e,Ae,Ve) , J^e,Hi^e, A„U \V,) 



m 



m 



m 



1 



> -{\^\{U,) + a\^i\{U2)) + j {p{xW + \h^l\^)dx + o{l). 
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Invoking T^ ni^Pe^ ^ej ^) > ^e ni^ei -^etU), we deduce that 
(4.8) 



lim inf ■ 



m 



>-{\f,\{U^) + a\ti\{U2)) + J {pix)\j\^ + \h-l\')dx. 



The left hand side of (|4.8p being independent from U, and Ui, U2 being arbitrary 
subsets of Si and 5*2 , we conclude that 

(4.9) limin f^-'^^^;^-'"^ >\f p{x)d\^^\ + f {p{x)\j\' + \h - l\^) dx . 

Now to conclude, we write 

^ = "^^^''^ + G' + ^^^''^) ' '^ = 'V + (/^ - /^m) • 

Upon substitution in the right hand side of l|4.9p and using in particular the re- 
markable identities 



curl ( J + TT-IT^^^M ] +h- hf_,=0 in Jl , /i^ - 1 = on c'fl , 



1 

p{x) 

we get the desired conclusion l|4.2p . 

Remark 4.3. Combining the upper and lower bounds of Propositions \3.1\ and 
Proposition \4-l\ then by uniqueness of the minimizer ^^ of E\ (see Section[^ below), 
it is evident that /i* = /i and h = h^^ . Here ^ and h are given in Proposition \4.1\ 
above. 



5. Minimization of the limiting energy 

As we explained in the introduction, by convexity and lower semi-continuity, the 
limiting energy (|1.5p admits a unique minimizer /^* which is expressed by means of 
the unique minimizer /i* of (|1.6p as follows, 



(5.1) 



^J■* 



-div 



1 



p{x) 



-V/i, 



Proceeding as in |12[ lllj. we may get an equivalent characterization of /i*. We 
write Hl{Q) for the space of Sobolev functions u such that u — 1 S HQ{n). 



Proposition 5.1. The minimizer u* of 



mm 

-div(^Vu)+ueM{n) 



A 



-div 



p{x) 



Vu +U+1 



is also the unique minimizer of the dual problem 



For instance, h^ 



u* + 1 minimizes the energy, 

mm / + 

feHl(n) \Jn P[x) 



2v\ dx. 



\f\'- 



p{x) 



dx , 



and satisfies —div I ■zrpi'^h^ j + /i* > 
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Proof. Let us define the lower semi-continuous and convex functional 

dx 



$(u) - ' ^(") 



n 2A 



-div ( \/u 



in the Hilbert space H = Hq{Q) endowed with the scalar product {f,g)H 
In pTxi'^f ■ ^-9 + fd- ^^t us compute its conjugate <!>*, i.e. 

$*(/)= sup {f,g)-^g). 

{g : *(3)<oo} 

Indeed, we have, 

$*(/) > sup / fijdx-— / p{x)\r,\dx- / fdx, 
riEL^ Jn ^-^ Jn Jn 

from which we deduce that 



**(/) = 



I +00 otherwise. 



By convex duality (see [HI Lemma 7.2]), 



min (\\u" 



, „2, + 2<i>(K)) = -min(j|/||2, + 2<i>*(-/)) , 

and minimizers coincide. That the measure //* = —div ( ^^V/i* 1 + /i* is positive 

is actually a consequence of the weak maximum principle, see [HI p. 131]. One may 
also follow step by step the proof given in [H]. □ 



Therefore, the limiting vorticity measure ^* is positive, and following [12], it can 



be expressed by means of the coincidence set wx = w\\J uP' 



(5.2) ^i={:,e^i : l-/,4^) = i-}, ^'={^^52:^^^ = ^}, 
as follows, 

(5.3) ^,, = (^i_^^l^^dx, 

where l^u;^da; denotes the Lebesgue measure restricted to w\. Furthermore, /i* (the 
minimizer of (|1.6p ) solves, 

-div (^V/i*) + /7, = inf^\wA 

/i* = 1 on dVt , 

and the regularity of /i» is complicated, as w\ may intersect the interior boundary 
dSi where /i, has gradient singularities (it is expected to satisfy a transmission 
condition, see the radial case below). But, we know that in the interior of S\ 
and 5*2, ft.* is locally C^'" for some exponent a G (0,1). Remark for instance 
that the measure /z* is no more uniform and may be discontinuous in light of 
/Lt* = (1 — ^) da; in w\ and /i* = (l — ^) dx in w\. 

For the sake of a better understanding of the sets wx, w\ and w\, we introduce 
the following critical constants (we emphasize their dependence on the material 
parameter a), 

(5.5) K{a)^ ./_, , .^ , Vie {1,2}, Ao(a)=min(Ai(a),A2(a)). 

2 max 



xes. 



p(a;) 
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Here we recall that ho is the solution of — div ( -ctV/iq ) + ho = in Q and ho = 1 
on dfl. The maximum principle gives that < ho < 1 in fl. We denote also by 
wo = w^U wl the set WAo(a) = w\^(a) U ^"Ao(a) introduced in dO]). 
The following proposition follows from a weak maximum principle [HI p. 131], see 
[T2] for a detailed proof (modulo necessary adjustments). 

Proposition 5.2. (1) w\ is increasing with respect to A and Ua>oWa = f^ ; 

(2) w\ and wf^ are disjoint; 

(3) If X < Ao(a) i/ien ft-* = /iq, M* = '^'^'^ ^''a = ; 

(4) If X= Ao(a) i/ien h^, — ho , fJ.* — and w\ — wq ; 

(5) // A > Ao(a) i/ien a^* 7^ ; 

(6) If X < Xi{a) for some i e {1,2}, then w\ = 0. 

Next for the sake of illustrating the above results, we give a detailed analysis of 
ho in the radial case. Let us take ft = D{0, 1) the unit disc in R^, Si = D{0,R) 
and ^2 = D{0, 1) \ D{0, R) where < i? < 1. In this case ho is radially symmetric, 
ho{x) = ho{\x\), and so it solves the following ODE, 



-h"{r)-r-^h'{r) + ho{r) =0 ii < r < R , 
-h"{r) - r-'^h'lr) + aho{r) =0 if i? < r < 1 , 
ho{R-) ^ ho[R+) , h'^{R^) = \h'^{R+) , 
h'o{Q) = fto(l) = l. 



(5-6) ^ r^^iu ^_ u^(T? \ unn \- i a' 



We look for a power series solution in the form, 

oo n 

ho{T)=^anr'' ifO<r<i?, /io(r) = ^ 5„(r - J?)" iiR<r<\. 

n=0 n=0 

The sequences (a„) and (&„) depend on a and R and all the terms are expressed as 
functions of the term ao, the two constants 

and the sequence (7„(a)) defined recursively by, 

7„3(a) = l, 7_2(a) = -^, 7„+i(a) = -— 7„(a) + a7„_i(a) , V n > -1 . 
Actually, we may verify that, 

a2„+i-0, a2„+2 = ^-^^-^^^^ VrieN, 

6o = aoa, 6i=aao/3, 6„+2 = 7 — — —7 (7„-2(a)&i + a7„_3(a)&o) V n > 2 , 

(n + 2)! 

and the term oq is expressed explicitly by. 



ao 



> ^(n + 2)! 



\ ?i=0 

provided that the sum in the r.h.s. is finite. For instance, this is the case when 
^ < i? < 1 and a ^ 0+. Actually, in the limit a — *■ 0+, it holds that, 

1 1-^0 1/,,, „, ii-a^ (-1)" (\ X\ 
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For i? e (5,1) chosen conveniently (close to ^), it is easy to verify that cq > 1 — q^^. 
Therefore, there exists oq G (0, 1) sufficiently small (depending on i? e (i, 1)) such 
that, for all a € (0, ao), we have, 

A2(a) < Ai(a), 

where Ai(a) and A2(a) are the critical constants introduced in l|5.5p . Coming back 
to the problem of vortex nucleation for the G-L energy (|1.3p . we get in light of 
Theorem 1 1.21 and Proposition 15.21 that for a wide range of applied magnetic fields, 

H = X\lne\il + oil)) (e ^ 0) , A G (A2(a), Ai(a)) , 

vortices exist and are pinned in ^2 . This result is in accordance with that obtained 
by the second author in [UIB], where nucleation of vortices near the critical magnetic 
field (in the case of the disc) is studied in details. 
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